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Summary 

Methodology has been developed for designing 
transonic airfoils and wings; it includes a technique 
that can account for static aeroelastic deflections. 
This procedure is capable of designing either super- 
critical or more conventional airfoil sections. Meth- 
ods for including viscous effects are also illustrated 
and are shown to give accurate results. 

The methodology developed is an interactive sys- 
tem containing three major parts. A design mod- 
ule has been developed that modifies airfoil sections 
to achieve a desired pressure distribution. This de- 
sign module works in conjunction with an aerody- 
namic analysis module, which for this study is a small 
perturbation transonic flow code. Additionally, an 
aeroelastic module is included that determines the 
wing deformation due to the calculated aerodynamic 
loads. Because of the modular nature of the method, 
it can be easily coupled with any aerodynamic anal- 
ysis code. 

Test cases are shown to demonstrate the viabil- 
ity of the method. The designs obtained generally 
match the target pressures and the airfoils used to 
generate them. For supercritical airfoil sections, the 
shock locations and the cusp regions were accurately 
reproduced by the method; however, some small dis- 
crepancies are noted near the leading edge. These 
differences are caused by the smoothing method cho- 
sen, and ways of eliminating them are discussed. 

Additionally, a case is shown that demonstrates 
use of the method to produce an entirely new design 
from an arbitrary target pressure distribution. This 
study found target pressure generation to be the most 
crucial step in the design process. 

Introduction 

A considerable amount of development work has 
been done in the area of transonic wing design. At 
present there is a wide spectrum of design codes avail- 
able, each with its own strengths and weaknesses. 
Each method is unique in some way, which makes 
grouping them difficult. The following paragraphs 
address some of the most often quoted classifications. 

One of the earliest design procedures, often clas- 
sified by itself, is the hodograph method by Bauer 
et al. (1972). This technique uses a variable transfor- 
mation to the hodograph plane which linearizes the 
partial differential equations for compressible poten- 
tial flow. 

Direct methods are another class of design ap- 
proach. These methods are characterized by a di- 
rect analysis on a wing or airfoil. A procedure is 
used to modify the shape, which is then analyzed 


directly to determine its new performance character- 
istics. Through several iterations of design estimates 
and performance predictions, a final design with the 
required performance is found. Often this process 
is automated, so that new design modifications are 
made by computer for each cycle. 

Shortly after the hodograph method was intro- 
duced, Barger and Brooks (1974) published a de- 
sign method based on streamline curvature. Their 
method of modification alters the airfoil to produce a 
pressure distribution that, with each iteration, more 
closely matches a specified target distribution. Davis 
(1979) has also published a direct method, which is 
an automated extension of the work by Barger and 
Brooks (1974). 

Optimization methods are also considered to be 
direct methods. They try to minimize a certain func- 
tion, such as drag, by modifying and analyzing a de- 
sign shape. Hicks, Murman, and Vanderplaats (1974) 
developed one of the first optimization methods for 
airfoil design based on Vanderplaats 1 (1973) opti- 
mization module, CONMIN. CONMIN continues to 
be widely used today. Kennelly (1983) has developed 
an airfoil design code using a quasi-Newtonian opti- 
mization method, QNMDIF. This method was shown 
to be more reliable and efficient than the conjugate- 
gradient method used in CONMIN. 

Inverse methods depend on a formulation of the 
design problem that starts with a desired (target) 
surface pressure distribution and then determines the 
airfoil shape that would have produced it. The tar- 
get pressures are used to determine velocities along 
the airfoil surface. These velocities are the bound- 
ary conditions for the flow field solution. After each 
iteration, the surface on which these boundary con- 
ditions are applied is changed in an attempt to re- 
duce the normal component of velocity. When the 
normal component of velocity is zero, the surface is 
the proper airfoil shape. Unfortunately, the airfoil 
surfaces derived may cross or leave an open trail- 
ing edge. One of the first transonic inverse methods 
was developed by Tranen (1974). Volpe and Mel- 
nik (1985) have formulated the first correctly posed 
inverse method for two-dimensional flow. Their tech- 
nique eliminates the problem of trailing-edge closure. 

Sobieczky et al. (1979) introduced the fictitious 
gas method. This formulation designs an airfoil that 
will yield shock-free flow at a given Mach number 
and angle of attack. His method alters the density 
of the gas in supersonic regions so that the flow 
equations remain elliptic. Then the supersonic region 
is recomputed from the potentials along the sonic 
line. The airfoil surface is altered to conform to a 
streamline through this supersonic zone. 
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More sophisticated analysis codes have brought 
many improvements to the early methods of airfoil 
design cited above. Many have been extended to 
include three-dimensional wings and viscous effects. 
Although many computational techniques now ex- 
ist for the automated design of transonic wings, a 
literature search shows that currently there are no 
published methods that consider the effects of acroe- 
lasticity. Deformation due to loads encountered in 
flight can have a significant effect on the performance 
of a wing by changing the local surface angle of at- 
tack, which causes a change in lift distribution. The 
procedure described in this document is the result 
of an effort to develop a methodology for transonic 
wing design that takes aeroelastic considerations into 
account. 

For this project, a direct method was developed 
that is similar to the work of Barger and Brooks 
(1974) and Davis (1979). The method is coupled 
with a small perturbation code to aerodynamically 
analyze complex configurations and with an aeroelas- 
tic module to include the effects of wing deflections 
on the performance of the configuration. 

A means of accounting for viscous effects is also 
included for use with the current method. Super- 
critical airfoil sections are especially subject to the 
decambering effects of a boundary layer, which can 
cause a significant lift loss at transonic speeds. 

Symbols 

b wingspan 

c local chord 

Cave average chord 

q, i = 1 , 5 smoothing equation coefficients 
C[ section lift coefficient 

Cp wing-body lift coefficient 

Cp pressure coefficient 

Cp local sonic pressure coefficient 

M Mach number 

R Reynolds number based on chord 

x streamwise coordinate 

y height coordinate 

a angle of attack 

rj spanwise location, 2y/b 

Description of Design Method 

Overview of Method 

A flowchart of the design procedure is shown in 
figure 1. The design method begins with the input of 


a starting shape and flow parameters to the aerody- 
namic analysis module. The shape information may 
include a single airfoil shape for a two-dimensional 
case, or a planform shape and airfoil sections for a 
wing. The aerodynamic analysis module is a tran- 
sonic small disturbance analysis code that employs 
an iterative solution procedure. Initially, several it- 
erations of the procedure are performed to obtain 
a rough estimate of the flow field. When the aero- 
dynamic analysis has completed a preset number of 
iterations through the solution, the method transfers 
control to the design module. This module redefines 
the shape of the airfoils based on the difference be- 
tween the current analysis and the target pressure 
distribution. This completes one design cycle. 

Each design cycle yields a modified airfoil or wing. 
The new shape is input to the aerodynamic analysis 
module where more iterations of the solution, usually 
about 25, are performed before returning to the 
design module for a new shape correction. The 
process is then repeated for a number of design cycles 
set by the user. 

Aeroelastic effects may be determined interac- 
tively, in conjunction with the design, or determined 
after the design is completed. When the interactive 
approach is chosen, the design method transfers con- 
trol to the aeroelastic module after every 20th itera- 
tion through the solution process of the aerodynamic 
analysis module (fig. 1(a)). The aeroelastic module 
determines the wing deflections and sends the de- 
formed wing shape back to the aerodynamic analysis 
module to be analyzed for 20 more iterations. This 
pattern allows at least one aeroelastic update to the 
analysis before each design shape is computed. Since 
the changes computed for each design cycle are rel- 
atively small, 20 iterations between aeroelastic cal- 
culations is more than sufficient. When aeroelastic 
effects are determined in the postprocessing mode, 
the aeroelastic module determines wing deflections 
once, after all other calculations have been completed 
(fig. 1(b)). 

Viscous effects can also be considered when using 
this design method. When the aerodynamic analy- 
sis module includes a boundary- layer calculation, it 
can be used interactively during the design process. 
Alternately, a boundary-layer thickness can be de- 
termined by postprocessing, and subtracted from the 
design. 

It is important to note that the success of the de- 
sign method is not dependent on the choice of aerody- 
namic analysis module. The wing- body- pod- pylon- 
winglet (WBPPW) transonic small disturbance code 
by Boppe (1987) was chosen for this exercise because 
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of its speed and ability to handle complex configu- 
rations, but the design method can be coupled with 
any aerodynamic analysis code. 

Design Algorithm 

The design algorithm employed was developed 
by Campbell and Smith (1987). It compares pres- 
sures computed by the analysis module with a 
target pressure distribution and computes a modi- 
fication to the airfoil section or sections. For three- 
dimensional cases, analysis pressures at the computa- 
tional stations are interpolated to the input locations 
and compared point by point with the input target 
pressures to determine the changes needed. The al- 
gorithm relates changes in surface pressure coefficient 
to changes in surface curvatures or slope, depending 
on local velocities. 

The algorithm originated with the observation 
from analyses of airfoils that changes in airfoil sur- 
face curvature generally gave proportional changes 
in the pressure coefficients. This relationship can be 
derived from the momentum equation in streamline 
curvature form, 

dq/q = -C(n) dn (1) 

where q is the local velocity, C is the streamline 
curvature, and n is the direction normal to the 
streamlines. Following the analysis of Barger and 
Brooks (1974), the distribution of curvature normal 
to the streamlines is assumed to be 

C(n) = C 0 e- kn (2) 

where A; is a constant and the subscript 0 denotes 
the value at the airfoil surface. For positive values 
of k , this expression yields the correct values of Cq 
and zero for the streamline curvature at the airfoil 
surface and in the far field, respectively. Substituting 
this expression into the momentum equation and 
integrating gives 

C 0 = k ln(<7o/?oo) (3) 

where the subscript oo indicates the free-stream 

value. If small changes in velocity are assumed the 
equation becomes 

A C 0 = k Ago/ <70 ( 4 ) 

Barger and Brooks (1974) suggest letting k be pro- 

portional to Co, which yields the following relation- 
ship between surface curvature and velocity: 

AC = A X C Aq/q (5) 


where A\ is a constant, and the zero subscripts have 
been dropped for simplicity. If small-disturbance 
assumptions are made (Ag — An, q — goo, Cp — 
-2u/goo, where u is the perturbation velocity in 
the free-stream direction and Cp is the pressure 
coefficient), equation (5) becomes 

AC - AC A C P (6) 

where A — -0.5Ai. In the above equation, as C 
approaches zero, a small change in curvature would 
result in a very large change in pressure coefficient. 
This does not occur in actual cases and probably 
results from the assumed normal curvatures distri- 
bution used in integrating the momentum equation. 
The curvature term has therefore been replaced with 
the following: 

(1 +C 2 ) B (7) 

where B is an input constant ranging from 0.0 to 0.5. 
This term approaches a value of 1.0 for small values 
of curvature, giving equation (6) the form 

AC = A A C P (8) 

For small values of surface slope where the curvature 
is approximately equal to the second derivative of 
the airfoil ordinate in the streamwise direction 
equation (8) is proportional to the expression given 
by Davis (1979), converted to incremental form, for 
subsonic flow past a wavy wall 

Ay" = k A C P (9) 

where k is a constant. For the larger values of cur- 
vature encountered in the nose region, the curvature 
term approaches C 2 ^. The value of B is reduced 
for cases with supercritical flow to help stabilize the 
solution. The final equation is therefore 

AC = A C P A{1 + C 2 ) B (10) 

with the sign of A reversing for the lower surface. 

While this algorithm works for both subsonic 
and transonic cases, it was found to converge rather 
slowly when strongly supercritical flow (local Mach 
number greater than about 1.15) was present. As 
pointed out by Davis (1979), the relationship of 
pressure coefficient and airfoil geometry would be 
expected to change character as the flow becomes 
locally supersonic. A hybrid algorithm similar in 
concept to Davis 5 was therefore developed based on 
supersonic thin airfoil theory. The equation for the 
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pressure coefficient is given in Liepmann and Roshko 
(1957) as 


C P = 2 y'J sjMl - 1 = k x y' (11) 

where y f is the surface slope. This can be rearranged 
to yield 

A y f = k AC P (12) 

for increments in Cp and y f . Differentiating this 
expression with respect to x gives 

Ay" = A;+(ACp) (13) 

Since the thin airfoil equation is strictly valid only for 
supersonic free-stream Mach numbers, an equivalent 
value for k must be determined. If the free-stream 
Mach number in equation (11) is replaced with the 
switching Mach number of 1.15, the value of k would 
be 0.28. In order to underrelax the changes in 
geometry, a value of 0.05 was used with good results. 
Several switching Mach numbers were also tried in 
various test cases. It was found that using the 
supersonic formula when the local Mach number just 
reached 1.0 adversely affected the design process. 
Apparently the supersonic character of the flow is 
not established until local Mach numbers of about 
1.15 1.20 are present. 

The hybrid algorithm is implemented by using 
the surface curvature relation for regions without 
strongly supercritical flow and a combination of the 
surface curvature and supersonic approaches in the 
higher Mach number regions. In order to modify 
the airfoil, changes in curvature are first converted 
to changes in y n using the formula 


Ay" = AC7 1 + 0/) 
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1.5 


(14) 


which assumes that the surface slope does not change 
appreciably. A change in y n at a single point can 
be made as shown in figure 2, To increase the 
magnitude of y n at point I without changing y n 
at the other locations, points I + 1 though N are 
sheared through a given angle. After all the points 
are modified as desired, the entire surface is rotated 
about the leading edge so that a given trailing- 
edge thickness is maintained (fig. 2(b)). This is 
important since crossed or open trailing edges could 
result if the leading-edge design is just slightly in 
error. This final rotation slightly changes the new 
values of curvature, but the changes are negligible 
as the solution approaches convergence. Because the 
method is applied point by point, a user can specify a 


particular region of an airfoil to be modified without 
disturbing the rest of the airfoil shape. 

The design algorithm is applied in the same man- 
ner whether the case is two- or three-dimensional. 
For a wing, the design modifications are determined 
and applied to each input section individually. The 
results of test cases have shown this approach to be 
satisfactory. 

Smoothing the airfoil section surfaces is found to 
be necessary to obtain a converged design solution 
for most cases. After the airfoil section is modified 
in each design cycle, the coordinates are smoothed 
by fitting them piecewise to two polynomials, one for 
the leading edge and a second for the remainder of 
the airfoil. Over most of the airfoil (aft of 20 percent 
chord), a third-order polynomial is used, 

y = c\ + C 2 X + c%x 2 + C4X 3 (15) 

with a least-squares fit. Sets of seven consecutive 
points along the airfoil surface are used to obtain the 
fit. Three points on either side of the center point 
determine the coefficients of the polynomial, and a 
new value is computed for the center point. Near the 
leading and trailing edges, the point to be adjusted 
may not fall in the center of the seven points. The 
leading- and trailing-edge points remain fixed, and 
the points near them are adjusted by fitting the curve 
through the nearest six points. The scheme is applied 
from the leading to the trailing edge so that all the 
points except the first and last are adjusted. 

In the leading-edge region, a third-order polyno- 
mial does not yield the best fit to the airfoil section; a 
square root term must be added to the polynomial to 
give the proper curvatures. This five-term equation 

y = c\ + C 2 X + C3X 2 + C4X 3 + C5X 0 ' 5 (16) 

can be unstable and cause a design to diverge when 
large changes to the airfoil section are required. 
Dropping the second- and third-order terms from 
equation (16) yields the three-term equation 

y = ci + C 2 X + C5X 05 (17) 

which improves the stability of the design method 
but can cause discrepancies between the design and 
target pressures. This effect is illustrated in test cases 
described later in this report. 

A procedure similar to the smoothing method is 
used to allow airfoil sections to rotate automatically 
during the design process. A rotation procedure is 
needed because the design algorithm is constrained 
to make its changes without moving the first and 
last points on each airfoil surface. As a consequence, 
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when angle-of-attack changes are needed in a design, 
the leading-edge curvature of one surface increases 
while the other decreases, resulting in a redistribu- 
tion of airfoil thickness, especially near the leading 
edge. Thus, the new camber or midthickness line will 
lie above or below the previous one, but will still be 
constrained to pass through the leading-edge point. 
By using the new camber line to extrapolate a curve 
through the leading-edge area, the airfoil coordinates 
can be adjusted in this area, and the section can be 
allowed to rotate (fig. 3(a)). 

In this procedure, the camber line coordinates 
are computed for all the input values of xjc in the 
leading-edge area, back to the sixth input point aft 
of xjc = 0.02. These last six camber line points are 
used to determine coefficients for a cubic polynomial 
curve fit through the points. The points ahead of 
xjc — 0.02, including the first point, are adjusted to 
lie on this curve. This adjustment is then applied 
to the corresponding airfoil coordinates to slightly 
reshape the leading edge (fig. 3(b)). Because the 
rotation changes are introduced to the procedure as 
small geometric changes, it is ensured that design 
changes to the surface curvature keep pace with 
the rotation changes and the design proceeds in a 
controlled fashion. 

Aerodynamic Analysis 

The wing-body-pod-pylon- winglet transonic small 
disturbance code (WBPPW) by Boppe (1987) was 
chosen as the aerodynamic analysis module for this 
study. It has the capability of modeling complex 
three-dimensional aircraft geometries with associ- 
ated pods, pylons, and winglets, as well as two- 
dimensional airfoils, at transonic speeds. The code 
has been used successfully for a wide range of con- 
figurations. A two-dimensional strip bound ary- layer 
approximation based on the method of Bradshaw and 
Ferriss (1971) is included. 

WBPPW uses two Cartesian grid systems in its 
solutions. The flow field is initially solved in a 
global crude grid system that encompasses the entire 
flow field. Embedded fine grids are used in areas 
of interest where flow field gradients may be large, 
such as near the wing and the fuselage. After the 
initial iterations have been computed on the crude 
grid, solutions are obtained on the fine grid, using the 
values at the crude grid points as Dirichlet boundary 
conditions on the outer boundaries of the fine grids. 
Then the crude grid values are recomputed, using the 
fine grid values as Neumann boundary conditions on 
the configuration surface. One pass through each 
of the two grid systems represents one aerodynamic 
iteration. A discussion of convergence considerations 
is included in the appendix. 


Structural Analysis 

The Transonic Aeroelastic Program System 
(TAPS) (Campbell, 1984) is used in conjunction with 
the aerodynamic analysis code and design module to 
complete the aeroelastic design method. It can be 
used interactively or in a postprocessing mode. The 
TAPS portion of the method determines the defor- 
mation of the wing shape under the current load. In 
the interactive mode that information is sent to the 
aerodynamic analysis code and updated after every 
20 crude/fine aerodynamic analysis iterations. In the 
postprocessing mode the deformation information is 
applied only to the final design. 

To use the TAPS module, a flexibility influence 
coefficient matrix must first be obtained for the wing 
structure. This matrix may be determined computa- 
tionally or experimentally. The coefficients produced 
are obtained by choosing “nodes” or points on the 
wing surface and then determining the vertical de- 
flection of each node when a unit normal force is ap- 
plied individually at each of the other nodes on the 
wing. This results in a square matrix containing, for 
every node, a deflection due to a unit force at each of 
the other nodes. In the TAPS module, the computed 
wing pressures are interpolated to the node locations 
and assumed to act uniformly over an area associated 
with each node. This area surrounds the node and is 
a quadrilateral whose edges lie half way between the 
current node and surrounding nodes. A force and a 
resulting deflection are then computed at each node 
location. These deflections are used to update the 
surface slopes, which serve as boundary conditions 
in the aerodynamic module. 

Results and Discussion 

Test cases have been performed to evaluate the 
accuracy of the design method. Several of them are 
discussed below. These cases demonstrate use of 
the procedure on supercritical airfoils and wings and 
show methods of including viscous and aeroelastic 
effects in the design process. Also, an applications 
case is presented, in which a new wing is designed 
from an arbitrary pressure distribution. 

All the test cases presented below were obtained 
using a standard procedure, with some variations for 
viscous or aeroelastic effects. The first step in prepar- 
ing a test case is to generate a set of target pres- 
sures. In all but the applications case, these pres- 
sures were generated by analyzing a known airfoil 
or wing with the WBPPW analysis code. Since the 
design method uses the same flow solver, differences 
between the design computation output and the orig- 
inal targets can be attributed to the design method 
rather than the accuracy of the aerodynamic analysis 
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procedures. After the target pressures were gener- 
ated, a NACA 0006 airfoil section was used as a start- 
ing airfoil and modified by the design method until it 
produced the target pressure distribution. Because 
large changes in the section shape were anticipated, 
the three-term equation (17) was used in the smooth- 
ing procedure. The results obtained are discussed 
below. For uniformity, all cases were allowed to iter- 
ate for 90 design cycles. Some of the cases converged 
more quickly, but 90 cycles have been found to be suf- 
ficient for almost all cases and are not prohibitively 
expensive. Solution convergence is discussed in more 
detail in the appendix. 

Supercritical Airfoil With Viscous Effects 

The first set of test cases illustrates use of the 
method to design a supercritical airfoil and shows 
methods for including viscous effects in the compu- 
tation. Viscosity can have a significant effect on 
the performance of a supercritical airfoil at tran- 
sonic speeds because of the decambering effect of the 
boundary layer. Figure 4 illustrates this by show- 
ing two analyses of the same airfoil, one with vis- 
cous effects included and one without. While there 
is little difference in the pressures corresponding to 
the forward region of the airfoil where the boundary 
layer is thin, further aft there is a large difference in 
shock location, as well as in the compression in the 
lower surface cusp region, and the viscous case ob- 
viously has much less lift. To ensure that an airfoil 
design will have the desired performance, one should 
account for viscous effects in the design procedure. 

Two procedures for including viscous effects are 
outlined in the following discussion. The first case 
determines viscous effects interactively, and the sec- 
ond uses a postprocessing method. 

Interactive boundary layer . The first ap- 
proach to viscous design is to utilize the interactive 
boundary-layer capability of the aerodynamic analy- 
sis code to couple the viscous and inviscid analyses 
during the design process. This allows the code to 
compute a boundary-layer thickness, apply it to the 
airfoil design, and have its effect included in the pres- 
sure calculations for each aerodynamic iteration. 

A section from the aeroelastic research wing 
ARW-2 (Seidel et al., 1985) was used in illustrat- 
ing this method. This wing is fully described below 
with the three-dimensional cases. The section was 
analyzed at Mach 0.77, an angle of attack of 1°, and 
a Reynolds number of 5 million to generate target 
pressures. Viscous effects were included in the anal- 
ysis. Then a design run was made, with the boundary 
layer being updated interactively during the design 
process. 


The results of this case are presented in figure 5. 
Figure 5(a) shows the target pressure distribution 
from the viscous analysis of the original supercriti- 
cal airfoil shape (target) as well as the pressure dis- 
tribution from the final cycle of the inviscid design 
run (final). Also included is the pressure distribu- 
tion computed from a viscous analysis of the designed 
shape (analysis). 

While the three distributions shown agree over 
most of the airfoil, the leading-edge pressures for the 
final and analysis computations are slightly differ- 
ent from the target pressures in that region. This 
discrepancy is a result of using the three-term equa- 
tion (17) in the smoothing procedure, as is illustrated 
in the section on smoothing considerations below. 

The resultant airfoil shape is shown in figure 5(b). 
It is compared with the original section used to gen- 
erate the targets. Note that the vertical scale is 
expanded to show detail. While there is a slight dis- 
crepancy in thickness, the overall shape and curva- 
ture appear to match very well. Too much curvature 
at some point on the leading edge of the lower sur- 
face causes the discrepancy seen on the lower surface. 
Also, a difference in the size of the trailing-edge thick- 
ness for the starting and target airfoils accounts for 
some of the discrepancy. 

It should be noted in this case and the following 
ones that the airfoil coordinates labeled “target” 
only define the shape used to generate the target 
pressures. It is expected that to reach the target 
pressures the design will need to replicate this shape; 
however, the design is in no way restricted to this 
shape. 

Postprocessing boundary layer . The ap- 

proach used in the postprocessing procedure is to 
compute the displacement thickness in advance and 
include it as part of the airfoil shape in an invis- 
cid design computation. This method requires sev- 
eral steps. First, a target pressure distribution must 
be determined that includes viscous effects. From 
this target pressure distribution, the boundary-layer 
thickness is evaluated. To ensure sufficient trailing- 
edge thickness, the boundary-layer thickness is added 
to the initial airfoil prior to the start of the design 
process. 

Once the initial shape has been determined, an 
inviscid design computation can be made. The air- 
foil generated by the design process will then include 
the thickness of the boundary layer and will need to 
be reduced by that amount to complete the design. 
This method requires less computer time than the 
previously described method and may be more ro- 
bust, as the sudden changes in thickness associated 
with boundary-layer development are eliminated. 
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The target pressures from the previous case were 
used for this design run. An estimate of the 
boundary-layer thickness was made from the tar- 
get pressures by using the method of Bradshaw and 
Ferriss (1971), which is part of the aerodynamic 
analysis code. This boundary-layer thickness was 
then added to the initial airfoil, a NACA 0006, to 
ensure sufficient thickness at the trailing edge. This 
initial shape was input to the design method, and 
an inviscid design computation was performed. The 
airfoil generated by the last design cycle was saved, 
and the previously computed boundary-layer thick- 
ness was subtracted from it. The resultant airfoil is 
the completed design shape. 

The results are shown in figure 6. Figure 6(a) 
shows the discrepancy in the leading-edge region due 
to smoothing as seen in the previous case. Addi- 
tionally, the analysis curve shows the shock to be 
slightly smeared, and the last 20 percent of the anal- 
ysis curve shows additional small deviations. Since 
these discrepancies are not present in the interactive 
case, they must be attributable to the bound ary- layer 
thickness calculation. The computed boundary layer 
has spikes of high negative curvature in the shock and 
trailing-edge areas. When these areas are subtracted 
from the smooth inviscid design, the resulting airfoil 
is not smooth. Several attempts were made at im- 
proving the correlation by smoothing the airfoil, but 
this was found to change the character of the airfoil 
in the process of eliminating the curvature spikes. 

Figure 6(b) shows the final design airfoil after the 
boundary-layer thickness has been subtracted from 
its upper and lower surfaces. This shape is not sub- 
stantially different from the interactive boundary- 
layer method design. Again the lower surface cur- 
vature has made the airfoil thinner than the original 
shape. Despite these small discrepancies, the results 
indicate that this is a viable method of accounting 
for viscous effects. 

This case required 166 seconds of cpu time on 
the NASA Langley Cray-2 computer for the inviscid 
design portion. For the previous case, 220 seconds 
were required. The time penalty is probably out- 
weighed by the simplicity and accuracy of the inter- 
active method. 

Interactive boundary layer with automatic 
twisting . The development of the automatic twist- 
ing procedure allowed its use in the third approach 
to this design case. Since the automatic twisting pro- 
cedure allows the leading-edge points to be changed, 
the curvatures in the leading-edge region are not as 
constrained as in the previous two cases. The inter- 
active boundary layer was also used for this case. 


In figure 7(a), the pressure distributions are 
shown to be well matched. There is a slight discrep- 
ancy near the trailing edge for the analysis run. Fig- 
ure 7(b) shows a very good comparison to the orig- 
inal airfoil. The thickness is well matched between 
the target and final airfoils. The only discrepancy 
is again in the trailing-edge area where the trailing- 
edge thickness is constrained to be slightly different 
for the two airfoils. 

Smoothing considerations . In each of the pre- 
vious cases, a discrepancy in the leading-edge region 
has been attributed to using the three-term equa- 
tion (17) in the smoothing method. To demonstrate 
this effect, the previous case was repeated using the 
five-term equation (16) in the leading-edge region. 
The results are shown in figure 8. The automatic 
twist procedure increases the stability of the five- 
term smoothing procedure; however, more design cy- 
cles are required to obtain a converged solution. This 
case required an additional 40 cycles. 

The pressure plot shows that the five-term 
smoothing method has eliminated the discrepancy in 
the leading-edge region seen in the three-term case. 
However, downstream the shock region match has 
been disturbed. Comparing figures 7(b) and 8(b) 
shows no significant difference between this airfoil 
and the one computed for the previous case. Again 
the trailing-edge thickness causes the aft portion of 
the airfoil to be slightly displaced from the targets. 

Plotting the surface curvatures of the original 
airfoil shows that the target section itself is not 
smooth. This could contribute to the discrepancies 
noted above, namely that the smooth design airfoil 
does not and can not match the original, rough 
airfoil. To test this theory, an additional design case 
was performed, this time smoothing the target airfoil 
shape before generating a new set of target pressures 
at similar conditions. Also for this case, the starting 
NACA 0006 airfoil surfaces were sheared slightly to 
increase the trailing-edge thickness to that of the 
design, eliminating any discrepancies due to that 
difference. Five-term smoothing is used to obtain the 
best match in the leading-edge region. The results of 
this case are shown in figure 9. 

The results of this case show the best agreement 
between design and target of all the cases in this 
section. The final design pressures pass through 
all the target pressures. The analysis only has a 
small discrepancy at the base of the shock. The 
airfoil upper surface matches all the targets; the 
lower surface has only a small discrepancy in the 
cusp region. This case demonstrates that the design 
method is capable of substantially reproducing a 
target shape and pressure distribution. 
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Supercritical Wing With Aeroelastic 

Effects 

After obtaining satisfactory results for the two- 
dimensional cases, this work was extended to three 
dimensions. The design algorithm was not modi- 
fied for three-dimensional applications; it was ap- 
plied in a streamwise, two-dimensional strip fashion. 
The calculated pressures, of course, reflect a three- 
dimensional flow field with sweep, root, and tip ef- 
fects included. The procedure was repeated for each 
input wing section, starting with the root section and 
proceeding out the span. 

Figure 10 illustrates the effect of aeroelasticity 
on a flow over a wing. The results shown are for a 
station near the tip of the ARW-2 wing. One pressure 
distribution was computed with aeroelastic effects 
considered; the other was computed using a rigid 
wing. For aft-swept configurations, static aeroelastic 
deflections typically result in negative twist angles, 
which decrease the lift obtained. 

In the following discussion, results from several 
test cases are shown. These illustrate two methods 
of accounting for aeroelastic effects, and the prob- 
lems of designing with no aeroelastic treatment. The 
first test case involves a design computation in which 
aeroelastic effects are considered interactively. The 
second case illustrates a design where aeroelastic ef- 
fects are not taken into account. The inconsistencies 
in twist angles result in an unacceptable solution. In 
the third case, this difficulty is resolved by allowing 
the sections to rotate during the design process. This 
leads to an alternate method of accounting for aeroe- 
lasticity, the postprocessing method. 

For the three-dimensional cases, the ARW-2 wing, 
a supercritical, aeroelastic wing (Seidel et al., 1985) 
was used. The ARW-2 is an aft-swept transport-type 
wing with an aspect ratio of 10.3, a taper ratio of 
0.4, and a leading-edge sweep of approximately 27°. 
A sketch of its planform is shown in figure 11. The 
10 streamwise lines across the planform in the figure 
indicate the stations where the design computations 
were performed. The fuselage was modeled in the 
computations as an infinite cylinder. 

Target pressures for all cases were generated by 
analyzing the wing in the WBPPW code at a Mach 
number of 0.80 and an angle of attack of 1.0° with 
the TAPS module operating so that aeroelastic ef- 
fects were included. A structural matrix was gen- 
erated for the ARW-2 wing using the Engineering 
Analysis Language (EAL) code (Whetstone, 1983), 
a finite-element, linear, structural analysis code. For 
the design work, the supercritical wing sections were 
replaced with NACA 0006 sections to start the pro- 
cedure; these sections were twisted to match the 
ARW-2 geometric twist distribution. The planform 


was held constant throughout the design process. 
Three-term smoothing was used. In order to con- 
centrate on aeroelastic effects, viscous effects were 
not considered in any of these computations. 

Case 1 — flexible design of ARW-2 wing . 
The first case includes the TAPS module in the de- 
sign run so that the aeroelastic deflections are cal- 
culated for each design iteration and are reflected in 
the modifications made to the wing. The results are 
shown in figures 12 and 13 for 3 of the 10 design sta- 
tions corresponding to the inboard, planform break, 
and tip locations. The curve marked “Final” shows 
the pressures computed during the final design cy- 
cle. This curve shows good agreement with the tar- 
get pressures at all stations. As shown previously, 
there is some discrepancy at the leading edge due 
to smoothing. The analysis curve shows pressures 
computed during an independent analysis of the de- 
signed wing. It shows the same agreement with the 
target pressures as was observed for the final design 
cycle. This demonstrates that the wing shape pro- 
duced does achieve the performance predicted by the 
design procedure. 

Figure 13 shows the designed airfoil sections. All 
the airfoils were slightly thinner than the original sec- 
tions used to generate the pressures, particularly the 
tip station. This is possibly a result of starting with 
the NACA 0006 airfoil, which has a thinner trailing- 
edge thickness than the ARW-2 sections. Because 
the first and last points on each surface are fixed and 
the designed surfaces must be sheared back to con- 
form to these points, the resulting airfoil is thinner 
than the original. Despite this difference, the surface 
slopes are closely matched. Even in the cusp region 
the contours are quite similar for the designs and the 
original sections. The analysis pressures seen in the 
previous figure demonstrate that these thinner air- 
foils are capable of producing the desired pressure 
distribution. 

Case 2 — rigid design of ARW-2 wing . The 
second case shows the results of simply designing 
to a pressure distribution when that distribution 
contains aeroelastic effects and there is no procedure 
in the design method to account for these effects. 
Figures 14 and 15 show the results at the same 
stations as the previous case. For the most part, the 
pressures match the target values closely, especially 
inboard, where less twisting would be expected to 
occur. However, at the tip there is a significant 
difference between target and design in the leading- 
edge region on the lower surface. The airfoils (fig. 15) 
follow the same pattern, matching closely inboard 
but at the tip showing a discrepancy. Since this is 
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essentially an angle-of-attack problem, matching the 
target pressures requires moving the stagnation point 
for the flow forward on the lower surface. This means 
decreasing the velocity and curvature at the lower 
surface leading edge and increasing the velocities and 
curvatures behind it. This results in the bulge seen 
on the lower surface. 

This design case demonstrates the necessity of a 
procedure for automatically changing the angle of at- 
tack of a section (twist). This led to the development 
of the automatic twisting procedure described previ- 
ously; the results of its application are shown in the 
following case. 

Case 3 — rigid design of ARW-2 wing with 
automatic twisting . The problem of case 2 was 
repeated for case 3. The same targets were used, and 
again no static aeroelastic deflections were calculated 
within the design computations (TAPS off). The 
same starting wing shape was used, but this time 
the automatic twist procedure described previously 
was employed. 

Figures 16 and 17 show the results of this run. 
The results in figure 16 are essentially the same as 
those for case 1, with the target pressures again 
closely matched except for the same small difference 
at the leading edge. The airfoils also show similar 
agreement after they are rotated to bring them in 
line with the target airfoils. The required rotation 
angles correspond to the deflection angles computed 
by TAPS in case 1. The differences between the 
rotation angles computed by the twist procedure 
and the deflection angles computed by TAPS were 
0.19° at the inboard station, 0.02° at the planform 
break station, and 0.20° at the tip. These results 
indicate that aeroelasticity need not be iteratively 
computed inside the design process. Twisting the 
wing during the design computations would be just 
as effective. This gives the user two modes for 
integrating aeroelastic twist into the design. 

One option is to use TAPS in the interactive mode 
to determine aeroelastic twist angles during the de- 
sign. This method is probably the most straight- 
forward, but it requires prior knowledge of the wing 
shape and structure for determining a flexibility ma- 
trix. This method would be most practical for mod- 
ifications to existing designs. 

The second option is to allow the sections to ro- 
tate as needed to reach the desired target pressures. 
This method would yield a shape and a twist distri- 
bution for the wing. The TAPS module could then 
be used in the postprocessing mode to compute the 
aeroelastic twist. The total rotation of each section 
would then be divided into aeroelastic and geomet- 
ric components. This technique would be the most 


versatile method and would be easiest to use on a 
completely new design. 

Application of Techniques 

The final test case is an application of the method 
to the task of designing an original wing from a set 
of arbitrarily defined target pressure distributions. 
The target pressures were developed to conform to 
certain design criteria, in contrast to previous cases 
where they were the result of an analysis of a known 
shape. Here, aeroelastic effects were taken into 
account using the postprocessing method described 
earlier. Viscous effects were determined interactively. 

In the initial step toward developing target pres- 
sure distributions for this case, a wing was fash- 
ioned using a natural laminar flow airfoil section at 
several stations along a slightly swept and tapered 
transport-type planform shape (fig. 18). The wing 
was mounted on an infinite cylinder fuselage. This 
configuration was then analyzed. The resulting pres- 
sure distributions were modified to conform to the de- 
sign criteria. These criteria included (1) maintaining 
a favorable pressure gradient over 60 percent of each 
airfoil, (2) maintaining a local Mach number below 
1.2 to delay boundary-layer separation and reduce 
wave drag, and (3) attempting to obtain an ellipti- 
cal span loading to reduce induced drag. In order to 
establish flow parameters, this design was assumed 
to be for a 10000-lb aircraft flying at 40 000 ft at a 
Mach number of 0.78. This yielded a Reynolds num- 
ber of 8.35 million, a dynamic pressure of 170 lb/ft 2 , 
and a design lift coefficient of 0.45. 

Once the target pressures were developed, they 
were used in a design computation. Again an 
NACA 0006 airfoil was used as the starting section. 
This section was combined with the previously de- 
scribed planform and body shape to yield a starting 
configuration. The automatic twist procedure was in 
use during this study. The wing design was begun 
with an angle of attack of 6.1° in order to obtain the 
necessary lift, but the twist angle of each individual 
section changed during the design process. 

Several iterations of target pressure modification 
and design work were required to obtain a satisfac- 
tory design shape. Target pressures had to be modi- 
fied to delay separation and to improve designs that 
were thicker outboard than inboard and airfoil sec- 
tions that were too thin or even had negative thick- 
ness near the trailing edge. 

Once a preliminary design was obtained, a flexi- 
bility matrix was computed for the wing, using the 
EAL method described earlier. For ease of com- 
putation, the wing was assumed to be solid, but 
built of a material flexible enough to give a real- 
istic tip deflection. Then further design work was 
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initiated, this time using the flexibility matrix com- 
puted for the intermediate design to allow interactive 
aeroelastic calculations. Once the design was final- 
ized, a new r flexibility matrix was computed and an 
aeroelastic analysis of the design was performed. 
Even though the design and analysis used different 
flexibility matrices, the twist angles computed by 
TAPS were nearly identical. The changes to the de- 
sign after the first flexibility matrix was computed 
were small and did not greatly affect the bending 
characteristics of the wing or the pressure differen- 
tials across the airfoil sections. 

Figure 19 shows the pressure distributions ob- 
tained in this design exercise. The targets are shown, 
along with the pressures from the final design cycle 
and the aeroelastic analysis of the final design shape. 
Figure 20 shows the final airfoil sections obtained. 
These results are shown at three stations along the 
wing: inboard, near the midspan, and near the tip 
(see fig. 18). 

The final and analysis plots show good agreement 
with the target pressures, matching the shock loca- 
tions and the pressure levels over most of the design 
at each location. Only some small discrepancies in 
the aft part of the distributions are noted. Design- 
ing realistic pressure distributions in this area was 
found to be a very difficult task. The designer must 
allow adequate airfoil thickness without causing sepa- 
ration. For the midspan section and the section near 
the tip, the upper surface leading-edge design and 
analysis pressures show a discrepancy with the tar- 
gets. This is due to smoothing. The results shown 
were obtained using three-term smoothing. Rerun- 
ning the case with five-term smoothing gave a much 
better match for the tip station, but for the section 
near the root gave a discrepancy in the leading-edge 
region of the lower surface. 

The span load distribution is shown in figure 21. 
The distribution obtained is not far from the design 
goal of elliptic. The analysis run results estimate the 
efficiency factor at 0.998. 

Some timing studies were done with this case. 
The work was done on a Cray Y-MP computer. The 
targets, analyses, and designs were run for 200 crude 
and 1000 crude/fine iterations. All runs included in- 
teractive boundary-layer calculations. For the initial 
design cases, where TAPS was not employed, 877 cpu 
seconds were required. When TAPS was included, in 
the later designs, 895 cpu seconds were used. The 
TAPS-on analyses used 721 cpu seconds. 


Figures 19 and 20 might be considered part of 
a series of working plots that a designer would use 
to obtain an acceptable wing design. While many 
important aerodynamic criteria were adhered to in 
developing this design, there may be other consid- 
erations that should be taken into account to incor- 
porate this into the total aircraft design. Balancing 
these considerations generally requires a number of 
iterations of modifications to the target pressures fol- 
lowed by design runs, so that developing an appropri- 
ate set of target pressure distributions becomes the 
goal of the designer. 

Concluding Remarks 

The effects of viscosity and aeroelastic deflections 
on wings at transonic speeds can be significant. A 
transonic design method has been developed for two- 
dimensional airfoils and three-dimensional wings; it 
includes techniques that can account for both vis- 
cous effects and aeroelastic deflections. Included in 
this work are suggestions for both interactive and 
postdesign techniques of accounting for these effects. 
This procedure is capable of designing either conven- 
tional or supercritical sections. 

Test cases are shown to demonstrate the viability 
of the method. Known airfoil sections and wings were 
analyzed to generate target pressure distributions, 
which were used as input for the design process. 
The designs obtained generally matched the target 
pressures and the airfoils used to generate them. 
Only small discrepancies are noted, usually near 
the leading edge where the smoothing procedure 
distorted the airfoil shape. The shock locations 
and the cusp regions were accurately reproduced by 
the method, even when viscous or aeroelastic effects 
were included. Interactive methods of accounting for 
viscous and aeroelastic effects were found to be the 
most accurate and easiest to use but required slightly 
more cpu time. 

Developing realistic target pressures for an arbi- 
trary design was found to be a difficult task. Many 
factors must be taken into consideration when devel- 
oping a target pressure distribution to ensure that it 
will yield a realistic design. The exercise is best ap- 
proached as an iterative procedure, using the results 
of design runs to guide modifications to the pressures. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
November 9, 1990 
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Appendix 

Convergence Considerations 

During development of the design method, it was 
noted that it was not unusual for a three-dimensional 
design case to require 200 crude and 2300 crude/fine 
grid iterations through the aerodynamic module be- 
fore a final design was obtained. This was many more 
iterations than were being used to determine the tar- 
get pressure distributions. A study was performed to 
determine how many iterations in the aerodynamic 
analysis module were required to adequately estab- 
lish the pressure distribution. In the following dis- 
cussion, iterations through the code are expressed as 
w a/6”. This reflects a initial crude grid iterations and 
b iterations using both the crude and fine grids. 

A section of the ARW-2 wing was used for this 
study. It was analyzed in the WBPPW code using in- 
creasingly larger numbers of iterations to determine 
the level at which the pressure distribution ceased 
to change. First, 100/80 iterations were used, con- 
forming to the recommendations of Boppe (1987). 
Next, this was increased to 200/300 and 200/1000. 
Finally, 200/2000 and 200/2400 were used to bracket 
the number of iterations typically used in a design 
case. 

For comparison purposes, a second wing was gen- 
erated by using the same planform as the ARW-2, 
but replacing its supercritical sections with a more 
conventional section, an NACA 64 a 012. This wing 
was also analyzed using the pattern of iterations 
listed above. Both three-dimensional cases were run 
at a Mach number of 0.8. The ARW-2 wing was run 
at an angle of attack of -1.0° and the convention 
section wing at 2.0°, to obtain a wing Cl near 0.55 
at the intermediate levels of iterations for both cases. 

Finally, a similar study was conducted for 2 two- 
dimensional airfoils, the NACA 64 a 012 section, and a 


supercritical section near the tip of the ARW-2. The 
two-dimensional cases were run at a Mach number of 
0.75 and a value of q near 0.55. 

The analysis code was found to require a higher 
number of iterations for convergence than was rec- 
ommended by Boppe (1987). Figure Al shows 
computed pressure distributions for the cases listed 
above. Figure Al(a) presents the results from the 
ARW-2 case, and figure A 1(b) shows the results ob- 
tained when the supercritical sections were replaced 
with conventional sections. For both cases, pressures 
near the tip (rj = 0.90) are shown. The figures clearly 
show that under these conditions the solution is not 
converged for the first two levels of iterations tested. 
There are significant differences between the magni- 
tudes of the pressures computed at each of the first 
three iteration levels. At the level of 200/1000 and 
above, however, the differences between the levels ap- 
pear to be small. For this reason, the target pressures 
for the three-dimensional example cases shown in the 
paper are computed with 200/1000 iterations. 

Figures A 1(c) and A 1(d) show convergence in 
the two-dimensional mode. Figure A 1(c) shows re- 
sults from an analysis of a supercritical section of 
the ARW-2 wing, and figure A 1(d) shows a two- 
dimensional analysis of the NACA 64 a 012 airfoil. 
The results here are even more dramatic. The dif- 
ferences in the solutions do not become small until 
the last two levels. 

The aerodynamic analysis code used in this 
method was found to require more iterations for con- 
vergence than was previously thought. This was 
taken into account, and target pressures were ob- 
tained at levels of convergence corresponding to the 
design procedure. For actual applications the large 
number of iterations required for a design will gener- 
ally ensure that the results are converged enough for 
engineering accuracy. 
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(b) 3-D conventional section; M = 0.8, a = 2.0°, C/, ~ 0.55. 
Figure Al. WBPPW convergence history study; n = 0.9. 




(c) 2-D supercritical section; M — 0.75, a = — 1.0°, q ~ 0.55. 



(d) 2-D conventional section; M = 0.75, a = 2.0°, q « 0.55. 


Figure Al. Concluded. 
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(a) With interactive mode. 



(b) With postprocessing mode. 

Figure 1. Automated predictor/corrector design method with aeroelastic effects. 
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(a) Pressure distribution. 



(b) Airfoil coordinates. 

Figure 6. Supercritical airfoil design case with postprocessing of viscous effects; M = 0.77, a = 1.0, 
R = 5.0 x 10 6 . 
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(b) Airfoil coordinates. 

Figure 7. Supercritical design case with interactive viscous effects and automatic twisting; M = 0.77, a = 1.0, 
R = 5.0 x 10 6 . 
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Figure 10. Effect of aeroelasticity on a section of a supercritical wing; M = 0.80; a 



Figure 11. ARW-2 wing geometry. 
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(c) Wingtip; 7? = 1.0. 
Figure 14. Concluded. 
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(b) Planform break; tj = 0.44. 
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Figure 16. Pressure distributions for the automatic twisting case; M — 0.80, a = 1.0. 
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(c) Wingtip; rj = 1.0. 
Figure 16. Concluded. 
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(b) Midspan; r) = 0.39. 
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(c) Outboard; T] = 0.95. 

Figure 20. Airfoil coordinates for the applications case; M = 0.78, a = 6.1; R = 8.35 x 10 6 . 
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Additionally, an aeroelastic module is included that determines the wing deformation due to the 
calculated aerodynamic loads. Because of the modular nature of the method, it can be easily 
coupled with any aerodynamic analysis code. 
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